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We present a measurement of time-dependent CP-violating asymmetries in decays of neutral B 
mesons to the final states D*^^ , using approximately 82 million BB events recorded by the BABAR 
experiment at the PEP-11 e + e~ storage ring. Events containing these decays are selected with a 
partial reconstruction technique, in which only the high-momentum ir from the B decay and the 
low-momentum tv^ from the D*^ decay are used. We measure the amplitude of the asymmetry to 
be -0.063 ± 0.024 (stat.) ± 0.014 (syst.) and compute bounds on | sin(2/3 + 7)|. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The Cabibbo-Kobayashi-Maskawa (CKM) quark- 
mixing matrix [1] gives an elegant explanation of CP 
violation and is under intense experimental investiga- 
tion aimed at overconstraining its parameters. A cru- 
cial part of this program is the measurement of the 
angle 7 = arg {-V ud V*JV cd V* b ) of the unitarity tri- 
angle related to the CKM matrix. The decay modes 
B° — > D* T h ± , where h = ir, p, on, have been pro- 
posed for use in measurements of sin(2/? + 7) [2], where 
(3 = arg {-V cd V* b /V td V; b )._ln the Standard Model the de- 
cays B° — > D* + tt~ and B° — » D* + n~ proceed through 
the b — ► Tied and b — > cud amplitudes A u and A c . The 
relative weak phase between the two amplitudes in the 
usual Wolfenstein convention for the CKM matrix is 7, 
which, when combined with B°B° mixing, yields a weak 
phase difference of 2/3 + 7 between the interfering ampli- 
tudes. 

The decay rate distribution for B — » D* ± -n T is 

e -|At|/r 

rfm = —ST- >< (!) 

[l ± S c sin(AmAt) =F ?7Ccos(AmAt)] , 

where t is the mean B° lifetime, Am is the B° — B a 
mixing frequency, and At is the difference between the 
time of the B — > Z)* ± 7r T (B roc ) decay and the decay of 
the other B (B tag ) in the event. The upper (lower) sign 
in Eq. 1 indicates the flavor of the B tag as a B° (B°), 
while r\ = +1 (—1) and ( = + (— ) for the B rec final state 
D*~ir + (D* + ir~). The parameters C and are 

Here S is the strong phase difference between A u and A c 
and r — \A U /A C \. Since A u is doubly CKM-suppressed 
with respect to A c , one expects r to be approximately 
2%. 

In this Letter we report a study of CP-violating asym- 
metries in B° — > D*^^ decays using the technique of 
partial reconstruction, which allows us to analyse a large 
sample of signal events. Additional information about 
the techniques used in this analysis appears in Ref. [3, 4]. 

The data used in this analysis were recorded with the 
BABAR detector at the PEP-II storage ring, and con- 
sist of 76.4 fb -1 collected on the T(4S) resonance (on- 
resonance sample), and 7.6 fb" 1 collected at an e+e~ 
center-of-mass (CM) energy approximately 40 MeV be- 
low the resonance peak (off-resonance sample). Samples 



of simulated Monte Carlo (MC) events with an equivalent 
luminosity 3 to 4 times larger than the data are analyzed 
through the same analysis chain. The BABAR detector is 
described in detail in Ref. [5]. 

In the partial reconstruction of a B° — > D* =F 7r ± candi- 
date (B rec ), only the hard (high- momentum) pion track 
TTh from the B decay and the soft (low-momentum) pion 
track 7r s from the decay D*~ — > D a ir~ are used. This 
eliminates the efficiency loss associated with the neutral 
D meson reconstruction. Applying kinematic constraints 
consistent with the signal decay mode, we calculate the 
four-momentum of the D, obtaining its flight direction to 
within a few degrees and its invariant mass (m m i ss ) [4]. 
Signal events peak in the m m i ss distribution at the nom- 
inal D° mass with an r.m.s. of 3 MeV/c 2 . 

In addition to B° — > Z?* =F 7r ± events, the selected event 
sample contains the following kinds of events: B° — > 
Z?* =F p ± ; BB background peaking in ™ m i S s, composed of 
pairs of tracks coming from the same B meson, with 
the 7r s originating from a charged D* decay, excluding 
B° — > D* T p ± decays; combinatoric B background, de- 
fined as all remaining BB background events; and con- 
tinuum e + e" — ► qq, where q represents a u, d, s, or c 
quark. We apply an event-shape cut and a D* helicity- 
angle cut to suppress combinatoric background. We re- 
ject TTh candidates that are identified as leptons or kaons. 
All candidates must satisfy 1.81 < m m ; ss < 1.88 GeV/c 2 . 
Multiple candidates are found in 5% of the events. In 
these instances, only the candidate with the m m i ss value 
closest to M D o is used. 

To perform this analysis, At and the flavor of the B tag 
must be determined. We measure At using At = (z roc — 
2tag)/(7/3c), where z roc (z tag ) is the decay position of the 
B r0 c (-Btag) along the beam axis (z) in the laboratory 
frame, and the e + e~ boost parameter 7/? is continuously 
calculated from the beam energies. To find z roc we fit 
the Wh track with a beam spot constraint in the plane 
perpendicular to the beams. We obtain z tag from a beam- 
spot-constrained vertex fit of all other tracks in the event, 
excluding all tracks within 1 rad of the D momentum in 
the CM frame. The At error UAt is calculated from the 
results of the z rec and z tag vertex fits. 

We tag the flavor of the B tag using lepton or kaon 
candidates. The lepton CM momentum is required to 
be greater than 1.1 GcV/c to suppress "cascade" lep- 
tons that originate from charm decays. If several flavor- 
tagging tracks are present in either the lepton or kaon 
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tagging category, the only track of that category used for 
tagging is the one with the largest value of Qt, the CM an- 
gle between the track momentum and the D momentum. 
The tagging lepton (kaon) must satisfy cos(9t < 0.75 
(cos#t < 0.50), to minimize the impact of tracks origi- 
nating from the D decay. If both a lepton and a kaon 
satisfy this requirement, the event is tagged with the lep- 
ton only. 

The analysis is carried out with a series of unbinned 
maximum-likelihood fits, performed simultaneously on 
the on- and off-resonance data samples and indepen- 
dently for the lepton-tagged and kaon-tagged events. The 
probability density function (PDF) depends on the vari- 
ables TO m iss, At, <TAti F, s t , and s m , where F is a Fisher 
discriminant formed from fifteen event-shape variables 
that provide discrimination against continuum events [4] , 
s t = 1 (—1) when the B tag is identified as a B° (B°), and 
s m = 1 (—1) for "unmixed" ("mixed") events. An event 
is labeled unmixed if the Wh is a 7r~(7r + ) and the -Btag is 
a B°(B ), and mixed otherwise. 

The PDF for on-resonance data is a sum over the PDFs 
of the different event types, V = "Y^nfiVi, where the 
index i = {D*ir, D* p, peak, comb, qq} indicates one of the 
event types described above, fi is the relative fraction of 
events of type i in the data sample, and Vi is the PDF for 
these events. The PDF for off-resonance data is V q q. The 
parameter values for Vi are different for each event type, 
except where indicated otherwise. Each Vi is a product 
of the PDFs M l (m miss ), T l (F), and T/ (At, a At , s t , sj. 

The ra m i ss PDF Mi for each event type i is the sum of 
a bifurcated Gaussian and an ARGUS function [4]. The 
Fisher PDF T% is a bifurcated Gaussian. The parameter 
values for J~d*-k, Fd*p, J~ pea k, and T com }, are identical. 

The At PDF is a convolution, 7? = 
/ dAt tT 7i(Ai tr , s t , s m ) Tli(At - Ai tr ,<7At), where % 
is the distribution of the true decay-time difference 
Attr and IZi is a three-Gaussian resolution function 
that accounts for detector resolution and effects such 
as systematic offsets in the measured positions of 
vertices [4]. 

The PDF T£,* 7r (At tr , s t , s m ) for signal events corre- 
sponds to Eq. 1 with 0(r 2 ) terms neglected, and with 
additional parameters that account for imperfect flavor 
tagging: 

7b. w = £ T (Att r ){a(l + a m /s) + (l-a)[(l-* t Aw) 
+s m (1 - 2uj) cos(AmAttr) 
-s t (1 - 2w) 5 s sin(AmAi tr ) ] } , (3) 

where the mistag rate ui is the probability to misidcntify 
the flavor of the B tag averaged over B° and B°, Auo is 
the B° mistag rate minus the B° mistag rate, a is the 
probability that the tagging track is a daughter of the 
signal D meson, k = 1 — 2p, where p is the probability 
that the daughter of the D results in a mixed flavor event, 
and £ = +(-) for s t s m = +!(-!). 



The -Btag may undergo a b — > ucd decay, and the kaon 
produced in the subsequent charm decay might be used 
for tagging. This effect is not described by Eq. 3. To 
account for it, we use a different parameterization [6] for 
kaon tags, in which the coefficient of the sin(AmAi tr ) 
term in Eq. 3 is 

- [(1 - ^) (* t a + * m c) + s t s m b(l - «»] , (4) 

where a = 2r sin(2/3 + 7) cos 6, b = 2r' sin(2/3 + 7) cos 5', 
and c = 2cos(2/3 + 7) (r sin S — r'sin<5'). Here r' (5') is 
the effective magnitude ratio (strong phase difference) 
between the b — > ucd and b — ► cud amplitudes in the tag- 
side decays. This parameterization is good to first order 
in r and r'. 

The PDF Tp* p for B° -> D*T p ± events is taken to 
be identical to T^,^, except for the CP parameters (S ± , 
a, b, and c), which are set to 0. The BB background 
PDFs T CO mb and T pea k have the same functional form as 
Eqs. 3 and 4, with independent parameter values. The 
parameters of T^ eak arc determined from a fit to the MC 
simulation sample. The CP parameters of all the BB 
backgrounds are set to and are later varied to evaluate 
systematic uncertainties. The PDF T q q for the continuum 
background is the sum of two components, one with a 
finite lifetime and one with zero lifetime. 

The measurement of the CP parameters proceeds in 
three steps: 

1. The parameters of M.i and the value of 
/-D*7r/(/r>*7r + Jd*p) are obtained from the MC simula- 
tion with the branching fractions B(B° — ► D*~ir + ) and 
B(B° — > D*~p + ) from Ref. [7]. Using these parame- 
ter values, we fit the data with Vi = -Mi(m miss ) Ti(F), 
to determine f comb , f D * p + f D *n, the parameters of 
A4qq, and the parameters of J-i for both continuum and 
BB events. This fit yields 6400± 130 (25160±320) signal 
events for the lepton- (kaon-) tagged sample. Projections 
of the fit results are shown in Fig. 1. The fit is repeated 
to determine the relative signal yields above and below 
the cut on cos 6>t, which determine the values of a and 
p. We find a = (1.0 ± 0.1)% ((5.6 ± 0.2)%) for lepton- 
(kaon-) tagged events. 

2. We fit the events in the sideband 1.81 < m m ; ss < 
1.84 GeV/c 2 to obtain the parameters of T comh and T' q -. 

3. Using the parameter values obtained in the previ- 
ous steps (except T^-), we fit the data in the signal region 
1.845 < TO m j ss < 1.880 GeV/c 2 , determining the param- 
eters of Tp t7r and TL. 

We use the MC samples to verify the entire analy- 
sis procedure, as well as the validity of using the same 
non-CP parameters in T^, p and T^,^ and of using the 
1~comb parameters obtained from the sideband in the sig- 
nal region. For lepton-tagged events, we find a bias of 
=F0.012 in S^, due to the assumption that events tagged 
with direct and cascade leptons are described by the same 
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mmiss(GeV/c 2 ) mmiss(GeV/c 2 ) 

FIG. 1: The m m i ss distributions for (a) lepton-tagged and 
(b) kaon-tagged events. The curves show, from bottom to 
top, the cumulative contributions of continuum, peaking BB, 
combinatoric BB, B° — > D* + p~ , and B° — » D* + tv~ events. 

resolution function. The results presented below are cor- 
rected for this bias. 

The signal-region fit (step 3) on the data determines 
the CP-parameters (S + , S~) for lepton tags and (a, b, c) 
for kaon tags which are found to be 

S + = -0.078 ±0.052 ±0.021, 

S- = -0.070 ±0.052 ±0.019, 

a = -0.054 ±0.032 ±0.017, 

b = -0.009 ±0.019 ±0.013, 

c = +0.005 ±0.031 ±0.017, (5) 

where the first error is statistical and the second is sys- 
tematic. The time-dependent, CP-violating asymmetry 
A C p = (N B o - N m )/(N B a + N m ) is shown in 

tag D Ug tag -°tag 

Fig. 2. In the absence of background and experimental 
effects, Acp = 2r sin(2/3 + 7) cos <5 sin(AmAi). 
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FIG. 2: The asymmetry Acp for (a) lepton- and (b) kaon- 
tagged events. The curves show the projection of the PDF 
from the unbinned fit. 

The systematic uncertainties on the CP-violation pa- 
rameters are summarized in Table I. They include 
(1) the statistical errors obtained from the fits of steps 
1 and 2; (2) uncertainties due to the unknown values 
of the CP-violation parameters in the background, the 
uncertainty in the ratio of branching fractions B(B° — > 



D*-ir)/B(B° -> D*~p), the modeling of T p ' eafc , and pos- 
sible biases introduced by the presence of background; 
(3) the uncertainty in the cascade lepton bias and pos- 
sible biases due to the r and Am parameters; (4) un- 
certainties in the measurement of the beam spot posi- 
tion, the detector z length scale, and detector alignment; 
and (5) the statistical error in the parameters determined 
from the MC sample. 

TABLE I: The systematic uncertainties on the CP- violation 
parameters. 



Source 


Error 


xlO 


- 3 ) 


n 




S + 


S~ 


a 


b 


c 


(1) Step 1 & 2 statistics 


1.7 


0.9 


1.0 


0.5 


0.6 


(2) Backgrounds 


12.1 


10.0 


13.7 


8.4 


14.2 


(3) Fit procedure 


6.6 


5.3 


5.2 


1.7 


0.8 


(4) Detector effects 


9.4 


7.3 


3.7 


9.1 


3.5 


(5) MC statistics 


12.8 


12.8 


8.0 


4.0 


9.0 


Total 


21 


19 


17 


13 


17 



Combining a and (S + + S )/2, accounting for corre- 
lated errors, we obtain 

2r sin(2/3 + 7) cos 6 = -0.063 ± 0.024 ± 0.014. (6) 

This measurement deviates from by 2.3 standard devi- 
ations. It can be used to provide bounds on | sin(2/3 + 
7)1 [8]. 

We use two methods for interpreting our results in 
terms of constraints on | sin(2/3 ±7)]. Both methods in- 
volve minimizing a x 2 function that is symmetric under 
the exchange sin(2/3 + 7) — > — sin(2/3 + 7), and applying 
the method of Ref. [9]. 

In the first method we make no assumption regarding 
the value of r. For different values of r we minimize the 
function \ 2 = fe=i ^jVjk 1 ^- x k, where Axj is the dif- 
ference between the result of our measurement and the 
theoretical value of S + , S~ , and a, and V is the mea- 
surement error matrix, which is nearly diagonal. The 
parameters determined by this fit are sin(2/3 + 7), which 
is limited to lie in the range [—1,1], and 5. We then 
generate many parameterized MC experiments with the 
same sensitivity as reported here for different values of 
sin(2/3 + 7) and with S — 0, which yields the most conser- 
vative limits. The fraction of these experiments in which 
X 2 (sin(2/3 + 7)) — Xmin i s smaller than in the data is in- 
terpreted as the confidence level (CL) of the lower limit 
on I sin(2/3 ± 7)|. The resulting lower limit is shown as a 
function of r in Fig. 3. This limit is always the more con- 
servative of the two possibilities implied by the ambiguity 
|sin(2/3 + 7)| <-> |cos<5|. 

The second method assumes that r can be estimated 
from the Cabibbo angle, the ratio of branching fractions 
B(B" D*tn-)/B(B° -> D*-n+) [10], and the ratio of 
decay constants fr>*j ' fo% [11], yielding r = 0.017i ;Q 07 . 
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FIG. 3: 95% CL lower limit on | sin(2/3 + 7 )| as a function of 
r. The solid curve corresponds to this analysis; the dashed 
curve includes the results of Ref. [12] for B° —* D*^n ± only. 

We attribute an additional non-Gaussian 30% relative 
error to the theoretical assumptions involved in obtaining 
this value. We minimize \ 2 = X 2 + A 2 (r), where A 2 (r) = 
for \r — ro|/ro < 0.3 and is an offset quadratic function 
outside this range, corresponding to a x 2 contribution 
with the uncertainties on ro given above. The parameters 
sin(2/3 + 7), 5, and r are determined in this fit. This 
method yields the limits | sin(2/3 + 7) > 0.87 at 68% CL 
and I sin(2/3 + 7 )| > 0.56 at 95% CL. 

Combining this measurement with the BABAR re- 
sults for fully reconstructed B° — > Z)* =F 7r ± and B° — > 
D^tt^ [12], we find, using the second method, | sin(2/3 + 
7)| > 0.87 at 68% CL and | sin(2/3 + 7) | > 0.58 at 
95% CL. The lower limit on | sin(2/3 + 7 )| obtained with 
the first method, including the results of Ref. [12] for 
B° — > D*^^ only, is shown in Fig. 3. 

To summarize, we have studied time-dependent CP- 
violating asymmetries in B° — » D^tt^ using partial re- 
construction. The amplitude of the CP-violating asym- 
metry that we measure differs from zero by 2.3 standard 
deviations. We interpret our results as an r-dependent 
lower limit on | sin(2/3 + 7 )|. 

We are grateful for the excellent luminosity and ma- 
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for the substantial dedicated effort from the computing 
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